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ABSTRACT 

Aims. We present the results of integral-field spectroscopic observations of the two disk galaxies NGC 3593 and NGC 4550 obtained 
with the Visible Multi Object Spectrograph at the Very Large Telescope. Both galaxies are known to host two counter-rotating stellar 
disks, with the ionized gas co-rotating with one of them. We measured in each galaxy the surface brightness, kinematics, mass surface 
density, and the stellar populations of the two stellar components as well as the distribution, kinematics, and metallicity of the ionized- 
gas component to constrain the formation scenario of these peculiar galaxies. 

Methods. We applied a novel spectroscopic decomposition technique to both galaxies, to disentangle at each position in the field 
of view the relative contribution of the two counter-rotating stellar and one ionized-gas components to the observed spectrum. We 
measured the kinematics and the line strengths of the Lick indices of the two counter-rotating stellar components. We modeled the 
data of each stellar component with single stellar population models that account for the of/Fe overabundance. 
Results. In both galaxies we successfully separated the main from the secondary stellar component that is less massive and rotates 
in the same direction of the ionized-gas component. The two stellar components have exponential surface-brightness profiles. In 
NGC 3593 they have diff'erent scale lengths, with the secondary one dominating the innermost 500 pc. In NGC 4550 they have the 
same scale lengths, but slightly diff'erent scale heights. In both galaxies, the two counter-rotating stellar components have diff'erent 
stellar populations. The secondary stellar disk is younger, more metal poor, and more of-enhanced than the main galaxy stellar disk. 
Such a diff'erence is stronger in NGC 3593 than in NGC 4550. 

Conclusions. Our findings rule out an internal origin of the secondary stellar component and favor a scenario where it formed from 
gas accreted on retrograde orbits from the environment fueling an in situ outside-in rapid star formation. The event occurred ^ 2 Gyr 
ago in NGC 3593 (1.6 + 0.8 Gyr after the formation of the main galaxy disk), and ^ 7 Gyr ago in NGC 4550, (less than 1 Gyr after 
the formation of the main galaxy disk). The formation through a binary galaxy merger cannot be ruled out, and a larger sample is 
required to statistically determine which is the most efficient mechanism to build counter-rotating stellar disks. 

Key words, galaxies: abundances - galaxies: individual: NGC 3593, NGC 4550 - galaxies: formation - galaxies: stellar content - 
galaxies: kinematics and dynamics. 



1. Introduction 

The general term of "counter-rotation" indicates the presence of 
two components of a galaxy that rotate along opposite direc- 
tions with respect to each other. This phenomenon has been ob- 
served in galaxies all over the Hubble sequence, from ellipticals 
to irregulars, including barred galaxies. Counter-rotating gala- 
xies are classed depending on the nature (stars vs. stars, stars vs. 
gas, gas vs. gas) and size (counter-rotating cores, rings, disks) 
of the counter-rotating c omponents (see Rubin 1994; Gallett^ 
119961: iBertola & Corshaj|[r999l for reviews). 

In this paper, we will focus on the class of disk galaxies 
with two counter-rotating stellar disks of comparable size. 
To da t e, the few known objects are: NG C 455 dRubin etaP 
[1992'; "Rix et al." ^1992"; "Em sellem et all [2001 . NGC 7217 
jMerrifield & Kuiiken 1 994), NGC 3593 jje rtola et al. 
[1996: iCorsini et al.i il998t. il999l: iGarcia-Burillo eTaLi |200Q) . 



Send offprint requests to: L. Coccato, e-mail: lcoccato@eso.org 

* Based on observations collected at the European Southern 
Observatory for the programme 087.B-0853A. 



NGC 4138 dJore et al.lll 996VHavnes et al.'l200Qh. and NGC 5719 
dVergani et al.ll2007l:ICoccato et al. 2011). 

The current paradigm that explains stellar counter-rotation 
is a retrog rade acquisition of e xternal gas and subsequent star 
formation dPizzella et al. I l200l . This scenario explains also the 
lower detection rate of counter-rotation in spirals with respect to 
SO's: a counter-rotating gaseous disk will be observed only if the 
mass of the newly supplied g as exceeds that of the p re-existing 
one dLovelace & Choulll996l: [Thakar & Rvdenlll996h . which is 
higher in spirals than in SO's. About 30 per cent of SO 's host a 
counter-rotating gaseous disk dBertola et al.| [r992: Kuii ken et al.l 
Il996t iBureau & Chung! 2006), although that less than 10 per 
cent of them have a s ignificant fraction of counter-rotating stars 
dKuijken et al.l Il996h . For comparison, less than 10 per cent 
of spiral galaxies are found to host a counter-rotating gaseous 
and/or stellar disk (Kannappan & Fabricant 2001; Pizzella et al.l 
12004 . 

Alternative scenarios for building counter-rotating stellar 
disks are also proposed. They include the accretion of already 
formed stars (and gas) through mergers and internal secular 
processes related to disk instabilities. Indeed, galaxy mergers 
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can produce galaxies with counter-rotating components, depend- 
ing on the orbital parameters of the encounter and the amount 
of gas involved in the process. The merger remnants range 
from kinematically decoupled core s or counter-rotating bulges 



in the case of minor mergers (e.^.. iBalcells & Gonzaleall99 



|Jesseit etal.ll2QQ7l: lEliche-Moral et al.1 |2Q1 ll: iBois et al.1 
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to large- scale ste llar disks in the case of binary majo r merg- 
ers (e.g. JPuerari & Pfenniger 2001; Cr ocker et al.l l2009l). An al- 
ternative to the e xternal-origin scenarios has been proposed by 
lEvans & CoUettl (1199 4) and it involves the dissolution of a triax- 
ial potential or a bar. In this process, the stars moving on box or- 
bits escape from the confining azimuthal potential well to move 
onto tube orbits. In non-rotating disks, there are as many box 
orbits with clockwise azimuthal motion as counter-clockwise. 
Thus, half box-orbit stars are scattered onto clockwise streaming 
tube orbits, half onto counter-clockwise. In this way, two iden- 
tical counter-rotating stellar disks can be built. Barred galaxies 
host q uasi-circular retrograde orbits (see Wozniak & Pfennige^ 
Il997h . Although t he origin of s t ellar counter-rotation ob served 
in barred galaxies (lBettonilll989l:lBettoni & Gallettalll997h is not 
necessarily external, acquired gas can be trapped on this family 
of retrograde orbits and then eventually form stars. 

The diff'erent formation mechanisms are expected to leave 
diff'erent signatures in the properties of the stellar populations of 
the counter-rotating component. In particular, the age is a key 
element to investigate the origin of counter-rotating stars. Gas 
acquisition followed by star formation predicts younger ages for 
the counter-rotating stellar component in all cases, and it al- 
lows for diff'erent metallicity and a-enhancement between the 
two disks, which are formed by diff'erent gas. Direct acquisition 
of stars through mergers also allows for diff'erent metallicity and 
Qf-enhancement, but the younger component will be determined 
by the diff'erence in age between the host galaxy and merged sys- 
tem. Roughly, one would expect younger counter-rotating stars 
in ~ 50 per cent of the cases. On the contrary, the internal origin 
predicts the same mass, chemical composition, and age for both 
the counter-rotating stellar components. 

A proper spectroscopic decomposition that separates the 
relative contribution of the counter-rotating stellar components 
to the observed galaxy spectrum is therefore needed. Such a 
decomposition will allow to disentangle the diff'erent forma- 
tion scenarios by properly measuring both the kinematics and 
the stellar population properties of the counter-rotating stel- 
lar disks. Unfortunately, the techniques that have been devel- 
oped over the years separat ed either the kinematics only (e.g., 
Kuiiken & Merrifieldl Il993l) or stellar populati ons only (e.g., 
Qd Fernand es et al.ll2005l:l Johnston et al.ll2012ah . 

Recently, we have presented a new technique of spectro- 
scopic decomposition which we have successfully applied to 
measure both the kinematics and stellar populatio n properties of 
the c ounter-rotating stellar disks of NGC 5719 (ICoccato et al.l 
I2OIII) . We found that the counter-rotating stellar component, 
which rotates in the same direction of the ionized gas, is younger, 
less rich in metals, more a-enhanced, and less luminous than the 
main galaxy disk. We therefore provided the crucial information 
that definitely confirms the gas-accretion scenario for NGC 5719 
previously suggested by Vergani e t al. ( 2007). 

Soon af t er, oth er groups developed similar codes. 
iKatkov et al.l (l2011bh appl ied an independent spectrosc opic 
decomposition technique (iKatkov & ChilingarianI 1201 lah to 
NGC 524, finding a stellar disk that counter-rotates with 
respect to the galaxy bulge. Both components are very old 
(15-20 Gyr); the disk component displays a steep metallicity 
gradient, while the bulge a more constant radial profile. Very 
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Fig. 1. Location of the VIMOS fields of view. The galaxy images 
are from the Digital Sky Survey. North is up and East left. The 
image scale is given on the bottom left corner in each panel. 



recently, 'Johnston et al.l (l2012bh analyzed a long- slit spectrum 
obtained along the major axis of NGC 4550 and found that the 
counter-rotating stellar component associated to the ionized gas 
is younger than the rest of the galaxy. 

In this paper, we apply our spectroscopic decomposition 
technique to two disk galaxies, NGC 3593 and NGC 4550, 
whic h are known to host two large cou nter-rotating stellar disks 
(Rub in et al.lll992l:lBertola et al.lll996b . We will: i) characterize 
the morphology, kinematics, mass surface density distribution, 
and stellar populations of both the counter-rotating components, 
ii) date the formation of the counter-rotating disks, and iii) probe 
their formation mechanism. We adopt a distance to NGC 3593 of 
7 Mpc ( Wiklind & Henkel 1992) and a distance to NGC 4550 of 
16 Mpc (Graham et al. 1999). The paper is organized as follows: 
in Section [21 we present the new integral-field spectroscopic ob- 
servations; in Section[3]we describe and apply our spectroscopic 
decomposition technique to disentangle the two counter-rotating 
components; in Section |4] we present the results; in Section [5] 
we discuss our results and compare them with the predictions of 
diff'erent formation scenarios; finally in Section Owe present a 
conclusive summary. 



2. Observations and data reduction 

The integral-field spectroscopic observations were carried out 
in service mode with the Very Large Telescope (VLT) at the 
European Southern Observatory (ESO) in Paranal during dark 
time between April 05 and June 02, 2011. The Unit Telescope 
3 was eq uipped with the Vis ible Multi Object Spectrograph 
rVIMOS. lLe Fevre et al.l[2Q03h in the Integral Field Unit (IFU) 
configuration. We used the 0.67 arcsec fiber"^ resolution and the 
HR blue grism, covering the spectral range 4150-6200 A with a 
reciprocal dispersion of 0.54 A pixel"^. The instrumental spec- 
tral resolution measured at 5200 A was 2.0 A (FWHM), equiva- 
lent to 115 km s"^. 

Observations were organized into diff'erent on-target expo- 
sures of 2950 s each, alternated to off'set exposures of sky fields 
of 280 s each for sky subtraction, as detailed in Table [T] For 
each galaxy we observed two overlapping fields-of-view aligned 
along the galaxy photometric major axis, at each side of the nu- 
cleus. The locations of the VIMOS fields of view are shown in 
Figure [T] 

Data reduction (bias subtraction, fiber identification and trac- 
ing, flat fielding, wavelength calibration, and correction for in- 
strument transmission) was performed using the VIMOS ESO 
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Table 1. Log of the spectroscopic observations. 



Target 


Field 


RA (J 2000.0) 


DhC (J 2000.0) 


OB 


Exp. time (+sky) 


Date 






[hh:mm:ss.s] 


[dd:mm:ss.s] 




[s] 


[yyyy-mm-dd] 
















NGC 3593 


A 


11:15:10.2 


12:56:53.7 


540380 


2850 (+279) 


2011-05-30 


NGC 3593 


A 


11:15:10.1 


12:56:55.1 


540382 


2850 (+280) 


2011-05-02 


NGC 3593 


A 


11:15:10.2 


12:56:55.8 


540383 


2850 (+280) 


2011-04-05 


NGC 3593 


A 


11:15:10.1 


12:56:54.4 


540384 


2850 (+279) 


2011-05-01 


NGC 3593 


B 


11:15:08.6 


12:56:55.1 


540385 


2850 (+279) 


2011-05-28 


NGC 3593 


B 


11:15:08.7 


12:56:55.7 


540386 


2850 (+279) 


2011-05-28 


NGC 3593 


B 


11:15:08.6 


12:56:54.4 


540387 


2850 (+279) 


2011-05-28 


NGC 4550 


A 


12:36:02.6 


12:05:32.5 


540388 


2850 (+279) 


2011-05-29 


NGC 4550 


A 


12:36:02.5 


12:05:31.9 


540390 


2850 (+279) 


2011-05-29 


NGC 4550 


A 


12:36:02.6 


12:05:33.2 


540391 


2850 (+279) 


2011-05-30 


NGC 4550 


B 


12:36:02.7 


12:05:10.5 


540392 


2850 (+279) 


2011-05-30 


NGC 4550 


B 


12:36:02.7 


12:05:09.9 


540393 


2850 (+279) 


2011-06-02 


NGC 4550 


B 


12:36:02.7 


12:05:11.2 


540394 


2850 (+279) 


2011-06-02 



Notes - Col. 1: target name; Col. 2: VIMOS field-of-view identifier, as given in Fig.[T] Cols. 3, 4: coordinates of the center of the field-of-view; 
Col. 5: observing block (OB) identifier; Col. 6: total exposure time on target (and on off'set position for sky background evaluation); and Col. 
7: date in which the OB was executed. 



pipeline version 2.6.3 under the EsoRex environmen fl plus 
some ad hoc IDL0 and IRAlQ scripts to correct fiber identifi- 
cation mismatch and remove cosmic rays. 

The diflTerent relative transmission of the VIMOS quadrants 
was corrected by comparing the intensity of the night- sky emis- 
sion lines. The sky oflTset observations were used to construct 
sky spectra. To compensate for the time variation of the rela- 
tive intensity of the night-sky emission lines, we compared the 
fluxes of the sky lines measured in the oflTset and on-target expo- 
sures. The corrected sky spectra were then subtracted from the 
corresponding on-target exposures. The (very) small amplitude 
of fringing patterns in the redu ced VIMOS spectra was corrected 
following the prescriptions by Lagerho lm et al.l (l2012h . 

Each exposure was organized in a data cube using the tabu- 
lated correspondence between each fiber and its position in the 
field of view. For each galaxy, the sky- subtracted data cubes 
were aligned using the bright galaxy nucleus as reference and 
co-added into a single data cube. Spectra from fibers map- 
ping adjacent regions in the sky were added together using 
the Voronoi binning method (ICappellari & CopinI l2Q03l) to in- 
crease the signal-to-noise ratio (S /N). Some of the spatial bins 
in NGC 3593 were modified to include only spectra from the 
regions associated with intense ionized-gas emission or the dust 
lanes crossing the galactic disk. We tested the robustness of our 
results using diflTerent binning schemes. 



Cappeflari & Emsefleml [20041) 



nalizea pixel titting coae (pr Ar, U 
developed in lCoccato et al.l (l201ll) . 

In each spatial bin, the code builds two synthetic templates 
(one for each stellar component) as linear combination of stellar 
spectra from the MILES library (ISanchez-Blazquez et al.ll2006h 
at FWHM = 2.54 A spectral resolution teeifi ori et aD l201ll) and 
convolves them with two Gaussian line-of-sight velocity distri- 
butions (LOSVDs) with diflferent kinematics. Spectra are nor- 
malized to their continuum level at 5 100 A, therefore the relative 
contribution of each component to the total spectrum is in terms 
of light. Gaussian functions are added to the convolved synthetic 
templates to account for ionized-gas emission lines (Hy, UJ3, 
[O III] ^^4959, 5007, and [N I] ^.15198, 5200) and fit simulta- 
neously to the observed galaxy spectra. Multiplicative Legendre 
polynomials are included to match the shape of the galaxy con- 
tinuum, and are set to be the same for the two synthetic tem- 
plates. The use of multiplicative polynomials also accounts for 
the eflTects of dust extinction and variations of the instrument 
transmission. 

The spectroscopic decomposition code returns the spectra 
of two best-fit synthetic stellar templates and ionized-gas emis- 
sions, along with the best-fitting parameters of luminosity frac- 
tion, velocity, and velocity dispersion. The line strength of the 
Lick indices of the two counter-rotating components will be then 
extracted from the two best-fit synthetic templates. 



3. Spectroscopic decomposition 

For both galaxies, we separated the contribution of the ionized- 
gas component and the two counter-rotating stellar disks to the 
observed spectrum using the novel implementation of the pe- 



^ The VIMOS data reduction pipeline and the ESO Recipe Execution 
Tool are available at http://www.eso.org/sci/software/pipelines/ 

^ Interactive Data Language is distributed by ITT Visual Information 
Solutions. It is available from http://www.ittvis.com/ 

^ Image Reduction and Analysis Facility is distributed by the 
National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 



3.1. Errors on measured parameters 

We carried on a set of Monte Carlo simulations on a set of 
artificial galaxy spectra to quantify the errors on line-of-sight 
velocity, luminosity fraction, and equivalent width as done in 
ICoccato et al.l (l201ll) . The simulations account for the character- 
istic velocity dispersion of our galaxies (< 90 km s"^), and the 
typical S /N ratios, wavelength coverage, and spectral resolution 
of our VIMOS observations. Errors in the measured quantities 
depend on the quantities themselves, and on their combinations. 

In the nuclear regions where the S/N is high (S /N > 50 
A"^), the errors are mainly aflTected by the small velocity dif- 
ference of the two stellar components (AV < 50 km s"^). In 
the outermost regions, where the velocity diflTerence between the 
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Fig. 2. Two-dimensional velocity fields of the main stellar component (left panel), secondary stellar component (central panel) and 
ionized-gas component (right panel) in NGC 3593. Scale and orientation are given in the left panel. Centers are in RA =11: 14:37.0 
(J2000.0), Dec. = +12:49:04 (J2000.0). 
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Fig.3. Same as Figured but for NGC 4550. Centers are in RA = 12:35:30.6 (J2000.0), Dec. = +12:13:15 (J2000.0). 



two stellar components is large (AV > 200 km s"^), errors are 
due to the low S /N (20 <S/N < 30 A"^). Errors on the line-of- 
sight velocity exceed 50 km s"^ for a luminosity fraction smaller 
than 20 per cent. Therefore, we set such a value as a lower limit 
for a component to be detected. In NGC 3593 errors in radial 
velocity do not exceed 40 km s"\ whereas in NGC 4550 they 
do not exceed 30 km s"^ Typical errors in the equivalent width 
are 0.6 A in NGC 3593 and 0.2 A in NGC 4550. 



4. Results of the spectroscopic decomposition 

We applied our spectroscopic decomposition technique to re- 
cover the kinematics, surface brightness, properties of the stellar 
populations, mass surface density of the counter-rotating stellar 
disks, and distribution and metallicity of the ionized-gas compo- 
nent of NGC 3593 and NGC 4550. 

In the following, we present our results. For each stellar com- 
ponent, we discarded from the analysis all the spatial bins where 
the contribution to the total luminosity is lower than either 20 per 



cent or 3(Tnoise where the noise standard deviation was computed 
from the residuals between measured and best-fitting spectra in 
the wavelength range 5000-5400 A. 



The data extend out to ~ 25'' along the photometric major 
axis of the two galaxies, which correspond to -0.8 kpc and ~2 
kpc for NGC 3593 and NGC 4550, respectively. The spatial sam- 
pling ranges from ~ 1.3 x 1.3 arcsec^ in the central regions up to 
-8x8 arcsec^ in the outermost bins. 



From now on, we will refer as "main" to the most massive 
stellar component, and "secondary" to the other. In the equa- 
tions, we will indicate them with the subscripts i and 2, respec- 
tively. We will show in Section l44l that in both galaxies the sec- 
ondary and less massive component is the one rotating in the 
same direction of the ionized gas. 
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Fig. 4. Maps (left and central panel) and radial profiles (right panels) of the surface brightness of the main and the secondary stellar 
components of NGC 3593 (upper panels) and NGC 4550 (lower panels). Scale and orientation of the maps are given in the left 
panels. The black, red, and blue solid lines correspond to the total surface brightness and to the surface brightness of the main and 
the secondary stellar components, respectively. The dashed red and dashed blue lines correspond to the best-fitting exponential disks 
to the surface brightness of the main and secondary stellar components, respectively. Dashed ellipses represent the boundaries of 
the elliptical annuli where the median surface brightnesses were computed. 



4.1. Stellar and ionized-gas kinematics 

The velocity fields of the main and secondary stellar compo- 
nents and ionized gas of NGC 3593 and NGC 4550 are shown 
in Figures [2] and O respectively. 

In the case of NGC 4550, the S/N was high enough to al- 
low the two LOSVDs to have diff'erent velocity dispersions. The 
S IN is lower in the case of NGC 3593, and therefore we decided 
to decrease the degree of freedom in the fit by setting the two 
LOSVDs to have the same velocity dispersion in the majority of 
the bins. This is a good approximation if we take into account 
the small diff'erence in velocity dispersion between the two stel- 
lar components measured by iBertola et al.l (Il996l) on long- slit 
spectrfl the spectral resolution, and the mean instrumental ve- 
locity sampling (^ 30 km s"^ pixel"^ at 5100 A) of the VIMOS 
observations. 

In both galaxies we did not observe a central peak in the 
stellar velocity dispersion; this indicates that the contribution of 
the bulge component is negligible. 



^ iBertola et al.l (Il996h measured cr ^ 60 km s~^ in the center and ^ 
70 km s~^ at 30'' where the main and the secondary stellar components 
dominate, respectively. 



Our results confirm the presence in both galaxies of a stellar 
disk and a ionized-gas disk that counter-rotate with respect to the 
main galaxy disk. 

NGC 3593 

The main stellar component reaches a line-of- sight rotation 
of Vi ^ 80 km s-i at 7? ^ 20'' (^ 700 pc) along a direction 
consistent with the galaxy photometric major ax is. This rota- 
tion is slightly higher than that measured by Be rtola et al.l (Il996l 
y ^ 70 km s"^) at the same distance. This is because in the pre- 
vious work the line of sight velocity distributions of the two stel- 
lar components were not fully resolved at 20''. The secondary 
stellar disk and ionized-gas components counter-rotate with re- 
spect to the main component, and their velocity is ¥2 - 100 
km s~\ consistent with that of the ionized-gas measured by 
IBertola et al] (11996 ^. 

The stellar velocity dispersion ranges between 30 < cr < 80 
km s"^ ; the velocity dispersion of the ionized gas is < 30 km s"^ . 

NGC 4550 

The main stellar disk rotates by Vi ^ 140 km s"^ at 7? ^ 25" 
(^ 2 kpc), while the secondary stellar disk and the ionized- 
gas component rotate slightly slower {¥2 - 110 km s"^). Our 
measurements are consistent with the long-slit absorption-line 
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Fig. 5. Equivalent width of the Lick indices in the main (red diamonds) and secondary (blue circles) stellar components in NGC 3593 
(upper panels) and NGC 4550 (lower panels). Predictions from single stellar population models by lThomas et al.l (l2Qllh are super- 
imposed. Crosses indicate mean error bars associated to the equivalent widths. 



and [Cj III] kinematics of iRix et al.l (Il992h and 'Johnston et al." 
(l2Q12bl) , and the ionized-gas kinema tics of [Sarzi et al. ( 2006) 
from HJS and [Oiii]. On the contrary. [Rubin et aP (ll992) mea- 
sured a larger rotation (V ~ 140 km s"^) for the Ha, [N II] 
and [S II] emission lines. Unfortunately, we do not have enough 
spectral coverage to investigate the difference in the kinematics 
of these emission lines. 

The velocity dispersions of the two stellar disks are very sim- 
ilar: the main component has 60 < cri < 80 km s"^ in the center, 
followed by a decline outside R > 14'' (1 kpc), where cri < 30 
km s"^ . The velocity dispersion of the secondary component in- 
stead ranges between 40 < 0-2 < 80 km s"^ The velocity dis- 
persion of the ionized gas is ^ 50 km s"^ . 

4.2. Stellar surface brightness 

We extracted the surface-brightness maps of the two stellar disks 
from their relative luminosity fraction (obtained in the spectral 
decomposition) and the observed total flux (obtained by collaps- 
ing the data-cube along the dispersion direction). 



The stellar surface-brightness maps and radial profiles are 
shown in Figure IH 

NGC 3593 

In NGC 3593, a preliminary fit with the IRAF task ELLIPSE 
(iJedrzejewskill 1 987h indicated that the two stellar disks have the 
same mean ellipticity and orientation on the sky within the er- 
rors. We therefore computed their median surface brightnesses 
along concentric ellipses, adopting the same ellipticity e = 0.55 
and position angle PA= 90°, both constant with radius. 

The surface-brightness profiles of the two stellar disks in 
NGC 3593 did not reveal any remarkable structure such as rings 
or spiral arms above the noise level (Fig. |4l). We parametrized 
their light distributions with two exponential disks ( FreemanI 
ri970) . which we found to have diff'erent scale lengt hs: = 
14"+2 " (480+70 pc) and /Z2 = 7"±1" (240+35 pc). BertolaetaLl 
(1996) reported similar scale length for the secondary compo- 
nent (h2 = 10", 340 pc), but larger value for the main compo- 
nent (hi = 40", 1.4 kpc). The discrepancy is due to the fact 
that we probe a diff'erent band and our data probe only the inner 
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2C, where the main stellar component is still faint. The differ- 
ent scale lengths of the counter-rotating stellar disks make the 
secondary stellar component to dominate the total galaxy light 
profile inside ^ 15", whereas the main stellar disk dominates 
outside. The transition radius between the light contribution of 
the two stel lar disks is consiste nt with the photometric decom- 
position bv lBertola et al.l (Il996l) . 

NGC 4550 

In NGC 4550, a preliminary fit with the IRAF task ELLIPSE 
indicates that the two stellar disks have similar orientation on the 
sky but diff'erent mean ellipticity. We therefore compute the me- 
dian surface brightness along concentric ellipses with the same 
PA= and diff'erent ellipticities (6i = 0.6, 62 = 0.5). 

The two stellar disks in NGC 4550 have the same scale 
length within errorbars {hi = 7fl ± Qf!4, i.e. 600 ± 30 pc, 
and h2 = 8V3 ± QfH, i.e. 650 + 50 pc) and the same luminos- 
ity, as found by Rix et al. (1992) from long-slit data. Their ap- 
parent flattening is difl'erent, meaning that they have difl'erent 
scale heights (Fig. [4]). Assuming the observed mean ellipticities 
(61 = 0.6 and 62 = 0.5, see Sect. 14.21) and a common inclination 
for the two stellar disks, we find that the minimum allowed incli- 
nation is / = 67°, which leads to intrinsic flattenings of ^0,1 =0 
and ^0,2 = 0.32. The intrinsic flattening qo was computed from 
^q(1 - cos^ /) = - cos^ /, for an oblate spheroid with measured 
flattening q =1-6 and inclination /. 

EmsellernetalJ (f2004) presented the mean stellar velocity 
field showing that the rotation is dominated by one compo- 
nent (corresponding to our main component) in a thin region 
along the major axis, while in a broad region away from the 
major axis the rotation is dominated by the other component 
(which rotates accordingly to the ionized gas). Their explana- 
tion is t hat the two componen ts must have very diff'erent scale 
heights (IC appellari et al. f l2007h . This is not in contrast with the 
mild diff'erence in the apparent flattening we found instead. In 
fact, the limitation of a single component fit is that it indicates 
only where one component is dominant; this leads to an overesti- 
mation of the flattening difl'erence. On the contrary, our spectro- 
scopic decomposition allows to measure both the two underlying 
spatial surface-brightness distributions, and to properly quantify 
the flattening. Alternatively, one can assume the same intrinsic 
flattening qo for both components. This would imply that the two 
stella r disks have difl'erent inclinations (lAfanasiev & Sil'chenkol 
l2002h . Extrapolating the typical value of SO/Sa galaxies (qo = 
0.18,|Guthrie 1992) to NGC 4550 (E7/S0), it results: h = 68° 
and 12 = 62° . 

4.3. Line-strength indices and stellar populations 

In each galaxy, we derived the properties of the stellar popu- 
lations of the two stellar counter-rotating components by mea- 
suring the line stren gth of the Li ck indices HyS, Mg b, Fe5270, 
and Fe5335 (Worthev^etaD ll994l) on the best fit synthetic tem- 
plates returned by the spectroscopic decomposition. To this 
aim, the spectra were set to rest-frame by adopting the mea- 
sured radial velocity and convolved with a Gaussian function to 
match the spectral resolution of the Lick system (FWHM = 8.4 
A. IWorthey & Qttavianilll997b . From the measured indices, we 
calculated the mean iron index (Fe) = (Fe5270 + Fe5335) /2 
("Gorgas et alj Il990l) and the combined magnesium-iron in- 
dex [MgFe]' = ^JMg b (0.82 • Fe5270 + 0.28 • Fe5335). The 
[MgFe]' index is almost independent fr om ^-enhancement and 
hence serves best as a metallicity tracer (iThomas et al ][2003|). 



Table 2. Luminosity-weighted values for the stellar population 
parameters of the steUar disks in NGC 3593 and NGC 4550. 





A^e 
[Gyr] 


rz/Hi 




NGC 3593 








Main: 


3.6 + 0.6 


-0.04 + 0.03 


0.09 + 0.02 


Secondary: 


2.0 + 0.5 


-0.15 + 0.07 


0.18 + 0.03 


NGC 4550 








Main: 


6.9 + 0.6 


-0.01 + 0.03 


0.20 + 0.02 


Secondary: 


6.5 + 0.5 


-0.13 + 0.04 


0.28 + 0.02 



Notes- [Z/H] and [of/Fe] are given in logarithm of solar units. Errors 
are computed as the standard deviation of the measurements divided by 
the square root of the number of spatial bins. 

In Figure [5] we compare the measured indices with the 
predictions of single-age stellar populations (SSP) model by 
iThomas et al.l (1201 1) that takes into account abundance ratios 
of a elements diff'erent from solar. Alternative line- strength in- 
dices as Hy^o (ICervantes & VazdekisI 12009 ) led to consistent re- 
sults. We decided to present the analysis based on the "classic" 
HyS sp ectral index for cons istency to our previous work on NGC 
5179 (Coccat o et al.]l2011h. 

We then fitted the I Thomas et al.l ( |201ll) models to the mea- 
sured indices to get two-dimensional maps of the luminosity- 
weighted mean stellar population age, metallicity ([Z/H]), and 
Qf-enhancement ([a/Fe]) of the two stellar components in the 
sample galaxies. They are shown in the left panels of Figure 
O and their luminosity-weighted values are given in Table [2l 
We also studied the radial dependence of these stellar popula- 
tion parameters, by computing their luminosity-weighted mean 
values along concentric ellipses using the same geometrical pa- 
rameters adopted in Sect. 14.21 They are shown in the right pan- 
els of Figure [6] Errorbars on the radial profiles are computed as 
<^ssp / Va^, where (Tssp is the standard deviation of the measure- 
ments within each ellipse, and N is the number of Voronoi bins 
within the ellipse. 

MILES stars belong to the Milky Way and this makes the 
stellar library trapped along the essentially one-parameter trend 
with a rise of [Mg/Fe] while decreasing [Fe/H] that translates 
into an anti-correlation between [a/Fe] and [Z/H] . This trend has 
a spread of ~ 0.5 dex in [Mg/Fe] and it is steeper in the metal- 
poor regime ( [Fe /H] < -0.4) while it is sh allower for stars with 
[Fe/H] > -0.4 (e.g.. iMilone et al.|[2QTl . Since we measured 
the Lick indices for linear combinations of MILES stars, the final 
values of [a/Fe] and [Z/H] we derived from SSP models could 
be not completely independent. 

The anti-correlation between [of/Fe] and [Z/H] is observed 
in the outer regions (R > 15'') of our galaxies where we mea- 
sured low values of metallicity (Fig.O, whereas it is less strong 
in the central regions (7? < 15'') where the metallicity is higher 
(0.15 < [Z/H] < -0.25). The ability of the fitting tool of re- 
covering the anti-correlation between (^-enhancement and metal- 
licity is an independent indication of its robustness. 

We used [a/Fe] as a proxy of the star-formation time scale A t 
following the approximation by (Thomas et al.l (l2005l Equation 
4). The main sources of errors on A^ are the uncertainties 
about SN la prog enitors and delay time distributions (e.g., 
iPritchet et al.ll2008h and dependence of our measurements of a- 
enhancement on metallicity. 

NGC 3593 
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NGC 4550 




Fig. 6. Stellar population parameters in NGC 3593 (upper panels) and NGC 4550 (lower panels). Left panels: maps of the age 
(top), metallicity (m iddle) and a-enhancement (bottom) of the two stellar counter-rotating components, as derived by fitting the 
iThomas et al.l (l2Qlll) models to the measured Lick indices. Dashed ellipses show the elliptical annuli used to bin the data and extract 
the radial profiles shown on the right panels. Orientation is the same as in Figures [2] and [3l Right panels: radial profiles of the 
total luminosity- weighted values within concentric ellipses of the age (top), metallicity (middle), and of-enhancement (bottom). Red 
and blue lines represent the main and secondary stellar components, respectively. The green line in the [Z/H] panels represents the 
metallicity of the ionized-gas component. 



The two stellar disks in NGC 3593 have diff'erent spectral 
properties. This diff'erence is visible in the diagnostic plots of 
Fig. [51 which compare the equivalent widths of Mg and 
Fe absorption lines. Unfortunately, the two-dimensional maps 
of the stellar population properties (age, [Z/H], [a^/Fe]) do not 
highlight this separation as clearly, because of the scatter in the 



measurements; the radial profiles in Fig. [6] serve better to this 
purpose. Although in some radial bins the two stellar compo- 
nents have consistent stellar population properties within the er- 
rorbars, there are clear systematic trends with radius. 

The age diff'erence between the two components is clearly 
visible for R < 15'' (~ 500 pc). In this radial range, the main 
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Fig. 7. Radial profiles of the mass-to-light ratio in the V band (left panels) and mass surface density profiles relative to the central 
value (right panels) for NGC 3593 (upper panels) and NGC 4550 (lower panels). Red and blue lines refer to the main and the 
secondary stellar components, respectively. The x-axis displays the semi-major axis of the ellipses shown in Fig.|4l 



component is the oldest (Agei(/? < 15 '0 = 3. 7+0. 6 Gyr) and the 
secondary component, associated with the gas, is the youngest 
(A^2(R < 15'0 = 1.4 ± 0.2 Gyr). 

The main component is systematically more metal rich 
([Z7H]i(5'' <R< 15^0 = -0.04 ± 0.04 dex) than the secondary 
component ([Z/UUS'' <R< 15^0 = -0.12 ± 0.03 dex). Both 
components have a metallicity peak in the central -150 pc: main 
component: [Z/H]i(7? < 5'') = +0.13 ± 0.04 dex, secondary 
component: [Z/UhiR < 5'') = +0.04 + 0.06 dex. The metallic- 
ity gradient of the secondary component is steeper than that of 
the main component. 

The main and secondary components are both characterized 
by an intermediate star-formation time scale (A^ ^2.5 Gyr) till 
15'' (^ 500 pc), whereas the outer parts of the secondary com- 
ponent are characterized by a much shorter time scale (A^ ^ 0.4 
Gyr). 

NGC 4550 

The spectral properties of the two stellar disks in NGC 4550 
are very similar, although there is an indication that the main 
component has slightly larger mean values of Mg b and Fe 
and smaller mean values of UJ3 than the secondary compo- 
nent (Fig. O. Our line- strength measurements of the Lick in- 
dices agree with the average values measured by lJohnston et al.l 



('201 2b|) for bo t h the co unter-rotating components. Nevertheless, 
John ston et al ] (I2012hh derived a mean age of the secondary 
component (2.5 Gyr), which is younger than ours (6.5 Gyr). This 
is due to the diff'erences in the adopted SSP models and per- 
formed analysis. Indeed, we used the models of Tho mas et al.l 
(2010) and computed the luminosity- weighted average of all 
the SSP parameter s we measured o n the tw o-dimensional field 
of v iew, wher eas [ Johnston et al.l (l2012bh fit the models of 
Vazdekis et to the arithmetic mean of the equivalent 

widths they measured in long- slit data. 

The total galaxy luminosity- weighted age within IRe - 15'.' 5 
(1.2 kpc) is Age= 6.9 + 0.4 Gyr. This is consistent with the lumi- 
nosity weighted age from a single component model (AgCKio = 
6.4^Qg Gyr, Kuntschner et al. 2010), although diff'erent authors 
claim ages larger than 10 Gyr (iAfanasiev & Sil'chenkol 120021) . 
The same is not true for the mean metallicity ([Z/H]= -0.04 + 
0.03 dex) and mean a/Fe abundance ratio ([cK/Fe]= 0.24 + 
0.0 3 dex) within l/?g , wher e luminosity weighted values found 
by iKuntschner et al.l (1201 Oh are more metal poor ([Z/H]kio = 
-0.25 + 0.04 dex) and slightly less a enhanced ([ck/FcIkio = 
0.15 + 0.05 dex). The diff'erence is probably due to the diff'erent 
iron lines used to probe the total metallicity and abundance ra- 
tio, and to the limited spectral resolution and wavelength range 
of the previous work. 
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Hg [ADU] ue [ADU] 

Fig. 8. Intensity map of the UJ3 emission line in NGC 3593 (left 
panel) and NGC 4550 (right panel). Fluxes are in arbitrary units. 
An asymmetric ring-like structure in NGC 3593 and a central 
asymmetric concentration in NGC 4550 are clearly visible. Scale 
and orientation are given. 
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Fig. 9. Map of the ionized-gas metallicity of NGC 3593 (left 
panel) and NGC 4550 (right panel). Dashed ellipses show the 
elliptical annuli used to extract the radial profiles of [Z/H] gas 
shown in Fig.[6l Scale and orientation are given. 



4.4. Mass surface density 

From the surface brightness and stellar population properties we 
can derive the mass surface density associated to each compo- 
nent. We used the stellar population models by Maraston ( 2005) 
to recover the stellar mass-to-light ratios from our stellar popu- 
lation measurements, and we then multiplied them for the light 
contribution measured from the spectra to get the mass. Since no 
flux calibration was performed, the mass profiles are relative to 
the central value. The radial profiles of stellar mass and mass-to- 
light ratios are shown in Figure H 

In both galaxies, the stellar component that rotates in the 
same direction of the gas is the less massive one. In NGC 3593 
this is true at all radii. In NGC 4550 this diff'erence is evident 
only in the inner 7'' (550 pc), whereas the mass distributions of 
the two stellar disks are very similar outside. 

4.5. Ionized-gas distribution 

The two-dimensional distributions of the ionized gas in NGC 
3593 and NGC 4550 derived from the fit to the 11/3 emission 
line are shown in Figure [8l These maps were produced using 
a specific Voronoi binning scheme with high spatial resolution, 
optimized to highlight the morphology of the HjS distribution. 

In NGC 3593, the map of the H/S emission displays an asym- 
metric ring-like structure with a peak-to-peak diameter of ~ 19'' 
(- 650 pc), which is more intense on the W side of the galaxy. 
The dimension and asymmetry of the Hp ring are ve ry similar to 
those of the HQf+[N II] ring found by Corsini et al. | (|T9 98), and 
to those of the CO ring found by Garcia-B uriUo et al.l (12000). 

In NGC 4550, the H/S distribution is patchy; in particular, 
there is a large double-peaked asymmetric structure in the central 
regions (either a ring or a disk) with a peak-to-peak diameter of 
~ 10" (~ 800 pc). The spatial extension of this H/S structure 
is similar to that of the regions where the mass density profile 
of the main disk is higher (see Section 14.41) and the secondary 
disk is younger (see Section 1431 right panel). The dimension of 
this inner I^ structu re is consistent with the region where CO 
emission is detected (ICrocker et al .112009b . 



The two-dimensional maps of the single stellar popula- 
tion parameters show a small systematic diff'erence in the total 
luminosity-weighted values of metallicity (see Tab.[2l), but not a 
significant diff'erence in age. This is due to the fact that the spec- 
tral HJS signatures are very similar for stellar ages greater than 4 
Gyr, and it is not easy to separate them with our accuracy (-0.2 
A). As in the case of NGC 3593, the stellar population values of 
the two stellar disks overlap within the errorbars in some radial 
bins, but systematic diff'erences are present. This is shown by the 
radial profiles in Fig.O the main component is more metal rich 
and less a-enhanced almost at all radii. 

The main diff'erence in age is observed in the central 5'' (400 
pc), where Agei(7? < 5'') = 7.7 ± 0.8 Gyr, and Age2(/? < 5'') = 
5.3 ± 0.7 Gyr. 

The main component has nearly solar metallicity at all radii 
([Z/H]i = -0.01 ± 0.03 dex), while the secondary compo- 
nent varies from nearly solar at the center ([Z/H]2(/? < 5'') = 
-0.03 ± 0.03 dex) down to subsolar outside 1.2 kpc ([Z/UhiR > 
15'') = -0.24 ± 0.05 dex). The metalficity gradient of the sec- 
ondary component is steeper than that of the main component. 

The derived mean stellar formation time scale of the main 
component is ^ 1 Gyr, and it is faster for the secondary com- 
ponent (A^ ^ 0.3 Gyr). 



4.6. Ionized-gas metallicity 

The gas-phase metallicity is one of the most important obser- 
vational diagnostics of the current evolutionary state of galaxies 
and, in the particular cases of NGC 3593 and NGC 4550, its 
measurements can give important information on the relation- 
ship between the ionized-gas and the two counter-rotating stellar 
components. 

We used the empirical relations between the intensity of the 
nebular lines and metallicity to estimate the gas-phase metal- 
licity. Since only the UJ3 and [Oiii] lines are present in the 
spec tral range we ad opted the R3 emission-lines ratio parame- 
ter (iLiang et al.ll2006h defined as: 

[O iii]^4959 + [O iii]^5007 

= h;^ • 

The galaxy spectra were not flux-calibrated. For this reason 
we followed the approach sugg ested by i Kobulnic ky & Phillip si 
( 2003) and recently adopted in Morelli e t al. (2012^ to replace 
the fluxes of the emission lines with their equivalent widths. This 
method does not need flux-calibrated spectra and has the further 
advantage of being insensitive to reddening. Thus, we measured 
the equivalent widths of the [Oiii]/l/14959, 5007, and BJ3 emis- 
sion lines in the best-fit spectrum. For each emission line, we 
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considered a central bandpass covering the feature of interest 
and two adjacent bandpasses, at the red and blue side, tracing 
the local continuum. The continuum level underlying the emis- 
sion line was estimated by interpolating a straight line in the 
continuum bandpasses of the spectra of the secondary compo- 
nent asso ciated with the g as. The used bandpasses are defined 
following Gonzalez (1993, Sect. 4.3) for HyS, [O 111] ^4959, as 
done in MorelH et al. (2012). 

The gas-phase metallicity measured for NGC 3593 and 
NGC 4550 was converted to [Z/H] following ISommariva et al.l 
(I2012h . The gas metallicity maps are shown in Figure |9l The 
luminosity-weighted mean values along concentric ellipses were 
obtained as in Sect. 14. 3 1 and are plotted in Fig.O 

NGC 3593 

The metallicity distribution of the gas shows an hint of a 
ring-like structure of similar dimension of the HJS ring detected 
in Fig. [HI The metallicities of gas and stars are consistent for 
R > 150 pc. On the contrary, in the nucleus of the galaxy, the 
gas metallicity remains constant at [Z/H]gas = -0.15 ± 0.09 dex, 
while the metallicity of the stars rises above 0.1 dex (see Sect. 
I4.3l and Fig. [6]). This discrepancy is probably due to the fact that 
in the very center of the galaxy the U/3 emission is almost absent 
(Fig. (8]); therefore in the center the star formation is less efficient 
today, and the present generation of stars in this region formed 
from a diff'erent enriched primordial gas. 

NGC 4550 

The gas metallicity distribution in NGC 4550 does not 
present particular structures in its two-dimensional map. Out 
to 20'' (1.5 kpc) the mean metallicity of the ionized gas is 
[Z/H]gas = 0.10 ± 0.07 dex, further out it peaks to 0.2 ± 0.2 dex. 
The gas-metallicity radial profile is in good agreement with the 
metallicity of the stellar components up to the very center, con- 
firming that the star-formation occurring in the gas disk is form- 
ing the secondary component. Contrary to NGC 3593, there are 
no diff'erences between gas and stellar metallicity in the center 
of the galaxy. This is due to the fact that in this galaxy the star- 
formation, as traced by the UJ3 emission, is occurring in a patchy 
structure all over the disk extension (Fig. [8]) and not mostly in 
the central structure. 



5. Discussion 

The morphology, kinematics, and properties of the stellar pop- 
ulations of NGC 3593 and NGC 4550 present both similarities 
and diff'erences, which help us in investigating their origin and 
assessing the most efficient formation mechanism. 

5.1. Origin of the large-scale counter-rotating stellar disks 

In each galaxy, the main stellar disk (i.e. the most massive) hosts 
a secondary stellar disk and a ionized-gas disk that counter-rotate 
with respect to it. In both galaxies, the stars in the secondary disk 
are on average younger, more metal poor, and more Of-enhanced 
than the stars in the main disk. The diff'erences in mass and 
age between the two counter-rotating stellar disks are more pro- 
nounced in NGC 3593 than in NGC 4550. In NGC 3593 the 
two stellar disks have the same intrinsic flattening but difl'erent 
scale lenghts, whereas in NGC 4550 the secondary component 
is slightly thicker than the main component, and they have the 
same scale length. Both galaxies host an asymmetric llj3 struc- 
ture in their central kpc of similar size: it is 650 pc wide with 



a ring-like shape in NGC 3593, and 800 pc wide in NGC 4550. 
These UJ3 structures match the regions where the secondary stel- 
lar disks are younger. 

Our findings support the scenario of a gas accretion followed 
by star formation as the origin for the stellar and gas counter- 
rotation in both galaxies. 

A pre-existing galaxy disk (that corresponds to the "main" 
component, as defined in Sect. ^ acquired gas from outside 
onto retrograde orbits. The amount of acquired gas must have 
been large enough to remove the pre-existing gas. The new gas 
first accreted in the outer regions of the galaxy, where formed 
rapidly new stars. The star-formation time- scale in the outer re- 
gions of the secondary stellar disk is indeed shorter than 500 
Myr in both galaxies, as inferred from the Of/Fe overabundance. 
Then, the star formation continued outside-in with a longer time 
scale, as inferred from the radial gradients in age and metallicity, 
which are more pronounced in the secondary component than 
in the main component. This process generated the secondary 
counter-rotating stellar components in both galaxies. These com- 
ponents are indeed younger than the corresponding main com- 
ponents and their angular momenta are aligned with that of the 
ionized-gas disks. 

The presence of molecular clouds, traced by the CO emis- 
sion, which co-rotate with both the secondary stellar disk and 
ionized-gas disk, is suggestive of star formation in the secon dary 
component (see Wiklind^ Henkel 19^2; Garcia-Burillo et al.l 
[2000 for NGC 3593 and lCrocker et al.l2Q09l for NGC 4550). The 
newborn stars have on average lower metallicity and, as a con- 
sequence of the anti-correlation between [Z/H] and [a/Fe], they 
are more Of-enhanced than the stars in the main stellar disk. This 
depends on the nature of the gas accreted by the two galaxies. 

The diff'erence between the luminosity-weighted ages of the 
two counter-rotating stellar components dates the gas accretion 
event. In NGC 3593 the accretion occurred between 2.0 and 3.6 
Gyr ago, i.e. 1.6 + 0.8 Gyr after the formation of the main stellar 
disk. In NGC 4550, our measurements do now allow a precise 
dating: it occurred ^ 7 Gyr ago, less than 1 Gyr after the forma- 
tion of the main stellar disk. 

Also the metal content of the ionized gas is consistent with 
the proposed gas-accretion scenario, although the fact that the 
metallicity of the gas component is similar to metallicity of both 
the stellar disks is intriguing. This may suggest that the acquired 
gas mixed with the pre-existing gas, responsible for the forma- 
tion of the main component. The final mixture of gas therefore 
depends on the relative enrichment and amount of the acquired 
and pre-existing gas, and on the additional enrichment subse- 
quent to the formation of the secondary stellar component. 

Although the metallicities of the ionized gas and the two 
stellar disks are similar on average, in the central 250 pc of 
NGC 3593 we measured a lower gas metallicity with respect to 
the stars (Sect. l4.6l) . This diff'erence between the two galaxies can 
be ascribed either to a diff'erent formation scenario or to a dif- 
ferent evolutionary stage in the formation process. In this latter 
hypothesis, the gas metallicity of NGC 4550 represents the final 
stage of the assembly process of the counter-rotating disk, which 
is still at work in NGC 3593: the stellar disk is forming outside- 
in, and the star formation will reach the center and then disap- 
pear, leaving stellar and ionized-gas components with approx- 
imately the same metallicity trend, as observed in NGC 4550. 
This interpretation is consistent with the diff'erent ages of the 
secondary stellar disks of the two galaxies. 

To date, the stellar populations of counter-rotating stellar 
disks of comparable sizes were successfull y disentangled only i n 
three disk galaxies. They are: NGC 5179 (ICoccato et al.ll201lh . 
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NGC 3593 and NGC 4550 in this work. In all of them, the stel- 
lar component co-rotating with the gas is younger, less massive 
and has different composition than the main stellar component, 
in agreement with the gas-accretion scenario. In NGC 5719, the 
gas-accretion event is traced by a stream of neutral hydrogen that 
connects the galaxy with its neighbor NGC 5713 (IVergani et al.l 
12007). Although we do not have the same evidence for NGC 
3593, its stellar populations are very similar to those of NGC 
5719; this makes us confident that these two galaxies share the 
same formation mechanism. The situation is not so clear for 
NGC 4550, as the diff'erences in the stellar population properties 
between the two stellar counter-rotating disks are not as evident 
as in NGC 3593. 

The scenario of internal origin of the secondary counter- 
rotating components is ruled out, as it would predict the same 
mass and stellar population for both components. The internal- 
origin scenario is ruled out even if we include a rejuvenation 
process for the secondary component (i.e. in the inner 850 pc of 
NGC 3593, and in the inner 400 pc of NGC 4550). In this alter- 
native case in fact, the secondary stellar component is expected 
to have higher metallicity (as the result of star formation from 
enriched gas) and lower cK-enhancement (as the result of the in- 
creased star-formation time scale, due to the coexistence of old 
and young stars) than the main stellar component in the regions 
where younger ages are measured. All these predictions are not 
consistent with our observations, and therefore we discard such 
a scenario. 

On the contrary, a binary galaxy merger cannot be com- 
pletely ruled out by our measurements as the formation mech- 
anism of counter-rotating stellar disks, given the limited galaxy 
sample and given that we have conclusive ev idence o f gas accre- 
tion only for one case. According to Cr ocker et al.l (1200 9), the 
merging scenario explains the fact that the gas rotates together 
with the thicker component in NGC 4550. On the other hand, we 
show that the diff'erence in thickness between the two counter- 
rotating stellar disks in NGC 4550 is not as pronounced as the 
one that can be guessed without a proper spectral decomposition. 
In addition, we found that the thicker stellar component is the 
youngest, and therefore it is the accreted component, whereas in 
the simulations the thicker component is the pre-existing stel- 
lar disk that >yas dy namically heated by the merging process 
(ICrocker et al.ll2009l) . 



5.2. The most efficient mechanism to form large-scale 
counter-rotating stellar disks 

We know that both binary mergers and episodes of gas accre- 
tion do occur in our Universe, therefore they are both viable 
mechanisms to create large-scale counter-rotating stellar disks. 
However, the relative fraction of counter-rotating components 
built by these processes is still unknown. 

As stated in the introduction, in the gas-accretion scenario 
the stellar component associated with the ionized gas is pre- 
dicted to be always younger than the main stellar component. 
In binary galaxy mergers instead, one would expect that the 
younger stellar component is associated with the ionized gas 
only in 50 per cent of the cases. Let's define M as the fraction of 
galaxies hosting two counter-rotating stellar disks of comparable 
sizes that were generated by a binary major merger (thus, 1 - M 
is the fraction produced by gas accretion), and E as the event of 
finding the stellar component co-rotating with the gas to be the 
youngest. 



The probability Pe to observe E is: 

Pe(M) = 1 - M + M/2 = 1 - M/2 (2) 

The probability 11^ to observe E in exactly N galaxies out 
of a sample of T counter-rotating galaxies is the first term of the 
Binomial distribution: 

n^(A^, r, Pe(M)) = (1 - Pe/-'^ (3) 

The value of M for which 11^ is maximum gives us the most 
probable value of the fraction of counter-rotating galaxies pro- 
duced by merger events. The condition: 

UE(N,T,PE(x))dx = 0.61 - I UE(N,T,PE(x))dx (4) 

M-M- Jo 

defines the Icr error bar on M. In our case (T = N = 3) we 
find that the fraction of counter-rotating galaxies built by binary 
mergers is M < 44 per cent. Therefore, measuring the diff'erence 
in age of the counter-rotating components in a larger sample 
of counter-rotating galaxies is fundamental to identify the most 
efficient mechanism. L arge spectrosc opic surveys like Atl as3D 
(iCappellarietal.ll201lb or CALIFA (ISanchez etal.ll2012h will 
help to identify other potential candidates by recognizing the 
kinematic signatures of stellar counter-rotating disks, like two 
symmetric peaks in the velocity dispersion map (iKrajnovic et al.l 
20111 

6. Summary 

We presented the VIMOS/VLT two-dimensional kinematics of 
the stars and the ionized-gas in the two disk galaxies NGC 3593 
and NGC 4550, which are known to host a stellar and a ionized 
gas disks that are counter-rotating with respect the main galaxy 
disk. 

We applie d the novel spectrosc opic decomposition technique 
introduced in Coccato et al.l (l201lh to separate and measure from 
two-dimensional spectroscopic data the morphologies, the kine- 
matics, the stellar populations, and the mass surface densities of 
both the two counter-rotating stellar disks in the sample galaxies 
for the first time. 

We found that in both galaxies the two stellar disks have dif- 
ferent stellar population properties, as measured by the equiva- 
lent width of the HyS, Mg and Fe absorption lines. The diff'erence 
is more evident in NGC 3593 than in NGC 4550. Both galaxies 
host a counter-rotating stellar disk, which rotates along the same 
direction of the ionized gas, and which is on average less mas- 
sive, younger, metal poorer, and more a enhanced than the main 
stellar galaxy disk. Our results support the scenario of external 
gas acquisition, followed by a subsequent outside-in star forma- 
tion as the origin of the observed counter-rotation. In NGC 3593, 
we are able to date the acquisition event that formed the counter- 
rotating stellar disk (1.6 + 0.8 Gyr after the formation of the main 
galaxy disk), while for NGC 4550 we can set only an upper limit 
(less than 1 Gyr after the formation of the main galaxy disk). 

Our results, combined with previous studies in the literature, 
rule out the internal scenario as the origin of counter-rotation in 
the studied galaxies. On the contrary, the merger scenario can- 
not be completely ruled out, given the small statistic available. 
Thus, a larger sample is necessary to identify the most efficient 
mechanism. 
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